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1. Introduction 



The object of our present study is to derive an expression for the optical conductivity as 
a convolution of the energy resolved joint density of states and an energy-frequency 
dependent transition rate. The need is to go beyond the usual reciprocal space 
based formulations and obtain an expression which we can immediately generalize for 
disordered systems. This would require labelling states by energy and the angular 
momentum labels (£, m) alone. Once we derive this expression we shall find a 
representation for the optical conductivity in the minimal basis set of the tight-binding 
linearized muffin-tin orbitals (TB-LMTO). The generalization to disordered systems will 
be carried out through the augmented space recursion (ASR) introduced by us earlier for 
the study of electronic properties of disordered systems [0-0- The ASR carries out the 
configuration averaging essential to the description of properties of disordered systems, 
going beyond the usual mean-field approaches and taking into account configuration 
fluctuations. The input into the ASR method includes the Hamiltonian parameters 
of the pure constituents, as the starting point of the local spin density approximation 
(LSDA) iterations for the alloy. It also includes the information about the transition 
rates of the pure constituents, expressed as functions of the initial and final state 
energies. The aim of this paper is to reformulate the reciprocal space representation 
of the transition rate and re-express it in the energy-frequency label representation for 
the pure constituents. Only when we are confident that this works, can we proceed 
with the full calculations for the disordered alloy. This communication is an attempt to 
verify our formulation for a series of alkaline earth titanates in the paraelectric phase, 
on which extensive theoretical and experimental data of optical properties are available 
for comparison. 

Perovskite structured titanate ferroelectric compounds is, to date, one of the most 
extensively investigated materials. They are extremely interesting from the viewpoint 
of solid state theoreticians because their structures are a lot simpler than that of any 
other ferroelectric material known, and, therefore, prove to be rather simple systems to 
study and better understand the ferroelectric phenomenon. The titanates are very easily 
prepared as polycrystalline ceramics, they are chemically and mechanically pretty stable 
and they exhibit para- to ferroelectric phase transition at or above room temperature. 
Continued interest in these compounds has led to a wide variety of theoretical and 
experimental work, specially on lattice vibrations. A less common approach has been 
those based on electronic structure calculations ||. Michel-Calendini and Mesnard 
have reported band structure of BaTiOs within a linear combination of atomic orbitals 
(LCAO) method with empirical off-diagonal integrals. The pioneering work on SrTiO^ 
was that of Kahn and Leyendecker || . This was followed by an Augmented Plane Wave 
(APW) calculation by Matheiss |10| and a self-consistent tight-binding calculation by 



Soules et al [11] . However, Battaye et al [p~2] have compared experimental valence-band 



spectra with these early theoretical predictions and have concluded that the agreement 
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was not satisfying. Pertosa and Michel-Calendini |13| carried out a modified tight- 
binding calculation on BaTiO^ and SrTiO% and compared their results with X-ray 
photoelectron spectra. These authors introduced inner orbital interactions. Perkins 
and Winter JT5) have carried out LCAO calculations on the band structure of SrTi0 3 . 
There have been several all-electron, full-potential linearized augmented plane waves 
(FP-LAPW) studies of the titanates in recent times [T7]-|T^|. In addition ultrasoft- 
pseudopotential, local density approximation (LDA) based studies on perovskites have 
been carried out by King-Smith and Vanderbilt |20| . In comparison, electronic structure 
calculations on CaTiO^ have been fewer. Ueda and coworkers [25|-[26| have used the 
first-principles tight-binding method to study CaTiO^ . 

We shall show that for all the three compounds the transition rate, defined by us, is 
strongly energy and frequency dependent, i.e. it depends upon the energy of both the 
initial and the final states. We shall compare the theoretical results with experiment. 



2. Methodology 

In recent years a number of methods have been proposed for calculating optical 
properties within the framework of the LMTO ||27||-|[32|| for both metals and semi- 
conductors. We shall present here a gauge-independent formalism, following the ideas 
of Hobbs et al ||32|| . Since our final aim is to use the augmented space recursion method 
(ASR) and study the optical properties of random systems, we shall modify the 

reciprocal space formulation and obtain an expression in which all states are labelled 
by their energy and the optical conductivity is expressed as a convolution of the energy 
resolved joint density of states and an energy-frequency dependent transition matrix. 
This formulation will then be directly generalized within the ASR. 

The Hamiltonian describing the effect of a radiation field on the electronic states of a 
solid is given by : 

H = E|^(p, + ^A(r i ,t)) 2 + V(r i ) + e$(r i ,t)| 

Here e is the magnitude of electronic charge, m e the electronic mass, c is the velocity 
of light and h is the Planck's constant. A(rj,t) and $(rj,t) are the vector and scalar 
potentials seen by the i-th electron because of the radiation field. There are N electrons 
labelled by i. The potential V^) experienced by the electrons is expressed as an 
effective independent electron approximation within the LDA of the density functional 
theory (DFT). For not too large external optical fields, neglecting terms of the order of 
0(|A| 2 ), the Hamiltonian reduces to : 

H = g{^P? + nr i ) + ij ( -A(r„ i )} (1, 
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Here jj = (e/m) pj is the current operator. We work in the Coulomb gauge where 
V • A(rj, t) — and $(rj, t) = 0, so that the electric field 

.9A(r t ,t) 
E(r^)- — 

In choosing the above equation we have ignored the response of the system. The local 
electric field is the external field due to the incident radiation as well as the internal 
field due to the polarization of the medium. Such local field corrections are important 
for insulators. We intend, as is customary, to introduce the local filed corrections as well 
as corrections due to the Coulomb hole in our final GW calculations, for which these 
single-particle picture will form the zeroth starting point. 

The Kubo formula then relates the linear current response to the radiation field : 

/oo 
dt> xAt-t')Mt') 
-oo 

The generalized susceptibility is given by : 

X^(r) = i6(r) (0o|[i M (r),i,(O)]|0 o ) 
where, r = t — t' and 0(r) is the Heaviside step function, 

0(r) = 



1 ifr > 
ifr < 



|0o) is the ground state of the unperturbed system, that is the solid in the absence of 
the radiation field. In the absence of the radiation field, there is no photocurrent, i.e. 
(0o I j mI^o) — 0. The fluctuation-dissipation theorem relates the imaginary part of the 
generalized susceptibility to the correlation function as follows : 



where, 



^» = \ (l - e-*>) S, v {u) (2) 



(5 = where ks is the Boltzmann constant and T the temperature 

k R T 



and 



/oo 
dt e lZT xA T ) z=^+t0 + 
-oo 



and, 



/oo 
dt e LZT (0 o |i M (r) >(O)|0 O ) z=cu+^0+ 
-oo 
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An expression for the correlation function, can be obtained via the Kubo-Greenwood 
expression, 

= ?EEE(^W/>-^/W^>^/-^-M ( 3 ) 

6 i f f 

We have assumed isotropy of the response so that the tensor is diagonal and we 
have defined S(u) as the direction averaged quantity | S^{uj). The are 
the occupied 'initial' single electronic states in the ground state while | {</>/}) are the 
unoccupied single electron 'final' excited states in the LDA description. 

The imaginary part of the dielectric function is related to the above : 

e 2 H = -L;S{u,) (4) 

We may obtain the real part of the dielectric function ex(u) from a Kramers-Kronig 
relationship. 

For crystalline semi-conductors the equation (|3[) may be rewritten as follows : 

= ^EE j BZ ^ K^-'kl J l*ik>l a 5(^(k)-^(k)-M (5) 

3 j' 

Here j and j' refer to band labels : j for the occupied valence bands and j' the unoccupied 
conduction bands at T=0 K. The k is the quantum label associated with the Bloch 
Theorem. For disordered materials, the Bloch Theorem fails and the expression (|5]) can 
no longer be used. Our first aim will be to obtain an alternative expression where the 
quantum states are directly labelled by energy and frequency, rather than by the 'band' 
and 'crystal momentum' indeces. For this, let us examine the following expressions : 



r a k 

occ unocc . 

J(E,u>) = EE / t-3 6 ( E ^) - E) 6(E r (k) -E-hu) (7) 

■/ J BZ o7T 
3 3 

In the equation (^), the right hand side picks up a factor of 1 whenever a quantum 
state, labelled by {k, j} falls in the range E, E+£E. The left-hand side, therefore, is the 
density of states arising from the bands labelled j. 

In the equation (0), the right-hand side picks up a factor of 1 whenever a quantum state 
in the filled bands labelled j falls in the range E, E+5E and simultaneously a quantum 
state in the unfilled bands labelled j' falls in the range E+cu, E+OJ+5E. The left-hand 
side is then the energy resolved joint density of states : 



J(E,u) = n v (E) n c (E + %uj) 



(8) 
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d 3 k 



T(E,u) 



r n k 

EE / T " ( k ) m<$ - E ) -E-fuj) 

■ / J BZ o7T 



3 f 



EE *(^-(k) - £) 5(%(k) -E-hu) 

■ i J BZ o7T 

(9) 

Where, 

T«'(k) = j M |*jk>| a 

The expression for S(u) from equation (||) then becomes, 

S(u) = (tt/3) / dE T(E,uo) J(E,uo) (10) 

Many earlier workers argued that the transition matrix element is weakly dependent on 
both E and u. They then assumed it to be constant T and obtained a simple expression 
for the correlation function : 

S (u) = (it/3) To J dE J{E,uo) (11) 

We shall investigate the validity of this approximation for the systems under study in 
this communication. Let us first get an expression for the equation fllQl) within the 
TB-LMTO formalism of : The basis of the LMTO starts from the minimal 

muffin-tin orbital basis set of a KKR formalism and then linearizes it by expanding 
around a 'nodal' energy point Efy . The wave-function is then expanded in this basis : 



**(r) = EEcia 



L a 



»+EE ^'(k) &(r) 



V a' 



where, L is the composite angular momentum index (£,m), j is the band index and a 
labels the atom in the unit cell. 

and, 



>vL\ 



t Y L (f) 
i l YAf) 



dE 



h%t,(k) = {Ct - E%) 5 LL ,6 aa , + J At S&(k) 
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C£ and A£ are TB-LMTO potential parameters and S^g^k) is the structure matrix. 
These terms are standard for the LMTO formulation and the reader is referred to the 



citation [21 1 for greater detail. The TB-LMTO secular equation provides the expansion 



coefficients c* La via : 



EE[^2'(k) + K - E*) 6 LL ,6 



V a' 



C L'a' 



(12) 



We may now immediately write an expression for the matrix element of the current 
operator as in equation (0) : 



<*i'k(r)| j |<& ik (r)> = <£, (Cl'WI j |# L (r)> • • • 

LL' a I 

... + {E*-E%) (Cl'WU |fe(r)>... 

... + {E^-E^) fc,(r)|j fc(r)>... 

... + {E*-E%) (E^-E^) fc,(r)| j fc(r)>} 

(13) 



We shall now obtain expressions for the right-hand terms by noting the following 



dr e 

J = e Tt = 

fffeW = ClM + Ci( r ) ( 14 ) 

We can write, 

r = (2tt/3) 1/2 r [(r w - F l5l ) i + i(Y 1 ^ 1 + Y Xjl ) j + 2 1 / 2 y 1 , k 
Using the above two equations we get, 

%v(t)*tH Cl Wrf 3 r = ^ f <,,(r) 0£(r) r 3 dr (15) 

Jo 

where, s Q is the atomic sphere radius of the a-th atom in the unit cell and is a 



combination of Gaunt coefficients [32 



I"/./.' = \l— 



In order to obtain / 0" L ,(r)* i/r 0° L (r) c? 3 r, we note that H = (h 2 /2m e )V 2 + V(r) so 
that using the Green's second identity we can obtain, 
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/fc(r)* H r C L (r) rf 3 r = 



^ LL' \ E^ £ , 



+ (h 2 /2m e )s 2 a <j>« e (s a ) <j,%,( 8a ) [PZ V - D% - 1)} 



(16) 



D" e is the logarithmic derivative of </>^(r) at r = s a and are obtained as parameters 
in the TB-LMTO routines. We define the following integrals : 



TO. 

1 nei 



Jo 



call, 



(R 2 /2me)4(^ - IKt> ~ 1) = ^ 
Then the matrix elements for the current operator becomes : 



(1) a 



LL' <f i 



ifvl = (€Ar)\ j, \€ L (r)) 



ifvl = (€Ar)\j,\€ L (r)) = 



T (4) a 
LL 1 ,ji 



= (€l^)\j, |fe(r)> = 







-i 




ft 




? /- 




-i 




n 






-e 


-i 




h 






-e 


-i 





1 LL> 



1 LL' 



E J$ £ + I pi — Vf, e <f>" e (s a ) (/>°gi(s a ) 
E Ku' + Jpe — Jeer ■ ■ ■ 

... — T>f, e <f>vi( s a) 4>'ii l { s a) 



1 LL> 



(17) 



where, £ = E° e - E° e , 

The transition term T-"'(k) has to be written in terms of the normalized wavefunction. 
The normalizing factor for the wavefunctions are obtained from : 



= J dh $* k (r) $ jk (r) 

= EE I4 k j 2 {J% + 2(E*-E: e )j% + (E*-E: e yj%} 



L a 



(18) 
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where, 



J La 

7 (2) 
J La 

T (3) 



Jo 

JO 



r 2 dr 



Using the secular equation ([12]), the expression for the transition term becomes, 



T jj '(k) 



-1/2 



E E E j'vJL + - *&) 42> • • • 



(19) 

The equation (10) provides an expression for the optical conductivity where both 
the transition matrix and the energy resolved joint density of states are expressed as 
functions of energy and frequency. As we shall show in a subsequent communication, 
that within the ASR formalism, this is the form in which the information about the 
constituents are input and the configuration averaged correlation function for the alloy 
may be expressed as : 



<S(cj)>= (tt/3) J dE T eff {E,u) <J(£,cj)> 



(20) 



where, 



< J(E, u) > = < n v (E) >< n c (E + uo) > 1 + A(E, to) < J(E, to) > 



and, 



T eff {E,u) = <T(£,cj)> +ST(E,(j,E{E, 



where, H(E,lj) is the self-energy due to disorder scattering and A(E,u) the 
corresponding vertex correction. The details of the derivation will be communicated 
in a subsequent paper |23 . 



3. Calculational Details and Results 

The primitive cell for the ideal perovskite structure ABC3 is illustrated in figure |I] 
(a). For the class of compounds we are interested in, the generic chemical formula is 
ABO3. A is a mono or divalent cation, B is a tetra or penta-valent metal. In the 
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(a) k z (b) 




k 

y 



Figure 1. (a) Cubic unit cell for a perovskite ABC3. (b) The Brilluoin zone for the 
cubic phase 

paraelectric phase there is full cubic symmetry. It can be thought of as lattice of corner 
sharing oxygen octahedra with inter-penetrating simple cubic lattices of A and B. The 
B cations sit in the centre of the octahedral O cage, while the A metal ions sit in the 
12-fold coordinated sites between the octahedra. In our case the body centre position is 
occupied by the Ti atom, the edges by alkaline earth atoms and the face centres by O 
atoms. The space group is 0\ and the corresponding Brilluoin zone is shown in figure [TJ 
(b). Both the A and B atoms are situated at sites with full cubic (Oh) point symmetry, 
while the O atoms have tetragonal (D±h) symmetry. 

The electronic configuration of the alkali earth atoms follows : 



Atom Deep Shallow Valence Unoccupied 

Core Core 



Ca 


Is 2 


2s 2 


2p 6 


3s 2 


3p 6 








As 2 




3d Ap 


Sr 


Is 2 


2s 2 


2p 6 


3s 2 


3p 6 


3d 10 


As 2 




5s 2 




Ad hp Af 


Ba 


Is 2 


2s 2 


2p 6 


3s 2 


3p 6 


3d 10 


As 2 Ap 6 Ad 10 5s 2 




6s 2 




5d 6p Af 


Ti 


Is 2 


2s 2 


2p 6 


3s 2 


3p 6 








3d 2 


As 2 


Ap 





Is 2 
















2s 2 


2p 4 


3s 3d 



Since we wish to take into account the shallow core states, to include the transitions 
from these to the conduction band at large enough optical frequencies, the energy range 
is about 40 eV (3 Ryd) and the single panel LMTO cannot be made to be accurate over 
this range, we have carried out a two panel calculation, with the lying in the lower 
energy range in one panel and in the upper energy range in the other. The minimal 
basis set used in the two panels are the following : 
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CaTi0 3 


Lower Panel 
Upper Panel 


Ca 4s Ap 3d, Ti 4s Ap 3d, 2s 2p (3d) 
Ca 4s Ap 3d, Ti 4s Ap 3d, 3s 2p (3d) 


SrTi0 3 


Lower Panel 
Upper Panel 


Sr 5s Ap Ad (Af), Ti 4s Ap 3d, 2s 2p (3d) 
Sr 5s hp Ad (Af), Ti 4s Ap 3d, 3s 2p (3d) 


BaTi0 3 


Lower Panel 
Upper Panel 


Ba 6s hp hd (Af), Ti 4s Ap 3d, 2s 2p (3d) 
Ba 6s 6p 3d (Af ), Ti 4s Ap 3d, 3s 2p (3d) 



The states in parentheses are unfilled states which have been downfolded in our 
calculations. The figure || shows the band structure of the three titanates. 



CaTi03 SrTi03 BaTi03 

25 H 1 1 1 h H 1 1 1 h H 1 1 




25 r Mxr Rr Mxr Ftr Mxr r 

Figure 2. Band structures of CaTiOs , SrTiOz and BaTiOz (No scissors operation 
has been carried out in these calculations) 

Let us first look at the leftmost figure for CaTiO% . The Ca 3p core level lies around 
-27 eV and not shown in this figure. The narrow band around -20 eV is the O 2s 
band. The nine valence bands just below the Fermi level are derived from hybridized 
Ti 4s and O 2p. An indirect band gap appears between the valence band top at the 
R point and the conduction band minimum at the V point. CaTiO% shows an indirect 
band gap (T-R) of 1.6 eV, in the absence of the scissors operation. The experimentally 
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reported indirect gap is 3.5 eV [ 25] . This discrepancy is characteristic of the Local 
density Approximation (LDA) upon which the TB-LMTO is based. In the conduction 
band region we have bands originating from (in ascending order of energy) Ti 3d t2 9 
triplet, a singlet arising from Ca 4s and a doublet from Ti 3d e g . Then comes the bands 
which originate from the Ca 3d e g doublet and the Ca 3d t 2g triplet. Finally we have 
the Ti Ap and Ca Ap based bands and finally the band based on Ca 4s. We note that 
for CaTiOz , Ca and Ti 3d based bands overlap and hybridize in the conduction region. 

For SrTiOz , in the lower panel, the Sr Ap level now sits almost atop the O 2s band, 
giving rise to a rather broad (as compared to CaTi0 3 ) s—p hybridized band just below 
-20 eV. The subsequent analysis of the bands is rather similar to CaTi0 3 . However, 
the band gap is now direct and ~ 1.4 eV. Earlier band structure calculations of Mo 
et al [[33|] and Kimura et al ]34| yield indirect band gaps of ~ 1.45 eV and ~ 1.79 eV 
respectively. The experimental direct band gap is around 3.2 eV. 

For BaTiOs , in the lower panel, the Ba 5p shallow core level now crosses and lies above 
the O 2s band. The band gap is direct and ~ 1.2 eV. The band structure is almost 
identical to the pseudopotential calculations of King-Smith and Vanderbilt p0| , whose 
band gap was also direct and ~ 1.8 eV. The experimental band gap turns out to be ~ 
3.2 eV |5|. 

For all three compounds, our calculations show a characteristic flatness of the lowest 
conduction band along T to X. This agrees with earlier works of Cardona [36], Matheiss 
ITJ , Harrison [TJ§, Wolfram and Ellialtioglu and King-Smith and Vanderbilt |20 



This observed flatness is related to certain unusual features in the density of states 
and optical conductivity, which appear to be characteristic of pseudo- two-dimensional 
systems. 

Figure ^| shows the densities of states for the three titanates. The densities of states 
reflect the detailed band structure we have described above. Earlier works on the 
density of states were based on different methods. Michel-Calendini and Mesnard 
|7], [8] and Pertosa and Michel-Calendini |13| have used parameterized tight-binding and 
adjusted LCAO based methods for BaTi0 3 . Their density of states is in agreement 
with ours, with the band gap difference characteristic of LSDA methods like the TB- 
LMTO. Similarly, Matheiss [|10j used the LCAO and Perkins and Winter [|I4[] used the 
extended Hiikel basis for SrTi0 3 . Again their band structures and densities of states 
are in agreement with ours, with the exception of the band gap. Even earlier works 
on SrTiOs by Zook and Castelman JTSI and Soules et al ]TB[| also show reasonable 
agreement. For CaTi0 3 we may compare our results with those of Ueda et al [25], 
which agree reasonably well with our results. 

Figure |] displays the transition rates T(E, u) for the three titanates, shown here as 
functions of the initial energy of the excited electron and the incident photon energy 
(frequency). It is clear from the figure that the transition rate is strongly dependent on 
the energy-frequency variables for all the three compounds. The usual assumption of a 
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CO 

O 



CO 

O 



CO 

o 





CaTi03 



SrTi03 





Figure 3. Densities of states for the three perovskite titanates 



transition matrix, weakly dependent on energy and frequency is certainly not valid in 
any of the three cases |3T|-|4"T|]. 

In figure |5| we compare the imaginary part of the dielectric function with the scaled joint- 
density of states/a; 2 . If the transition rate were independent of energy and frequency, 
they should be the same. The behaviour of the two are similar, but the relative weights 
of the structures across the frequency range are clear indications of the energy-frequency 
dependence of the transition rates. 



CaTiO? 



The effect of energy-frequency dependence of the transition rate has a large effect for 
the optical properties of CaTiO^ . If we compare the results reported in with 
our figure |^ we note that although in the earlier work the joint density of states does 
reproduce the peaks at lower frequencies, the relative heights are not replicated. Looking 
at the lower panel of the figure, we may assign the peak at 4.5 eV to the transitions : 
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Figure 4. Transition rates for the three perovskite titanates shown as functions of 
initial energies and incident photon energies. BaTiO^ (top), SrTiO^ (left-bottom) 
and CaTiO^ (right-bottom) 

O 2p — > Ti 3d-t2 g at the R and X-points. The next and highest peak arises because of 
the nearby two unresolved peaks due to the transitions : O 2p — > Ca 4s at the M-point 
(at 6.3 eV) and O 2p —>■ Ca 3d-e g at the M-point (at 6.8 eV). The third peak at 7.5 
eV may be assigned to the transitions : O 2p — > Ti 3d-t2 9 at the R-point and O 2p 
^Ca 3d-e g at the X-point. At higher frequencies, the theory does not tally well with 
experiment. The peak at 10 eV is not reproduced except as a shoulder. The real part of 
the dielectric function ei(u) is obtained by a Kramers-Kronig transformation from the 
imaginary part. This is shown in the top panel of figure [| The discrepancies at high 
frequencies could be due to the fact that in the scissors type approach we have provided 
a rigid shift to the conduction bands. In a full fledged many-body GW technique, which 
is our ultimate aim to produce, the shift, due to the self-energy will turn out to be 
energy (frequency) dependent. 

SrTi0 3 : 

Let us now examine the figure [7]. As in the case of CaTiO^ , here too, the effect 
of energy-frequency dependence of the transition rate reproduces the correct relative 




Figure 5. Comparison of the calculated imaginary part of the dielectric function (€2) 
(left) for the three perovskite titanates with the same function calculated considering 
transition rate is constant (right). 



heights of the peaks in e 2 (uj). The shoulder peak at around 4 eV may be attributed to 
the transitions : O 2p at — > Ti 3d — t 2g and O 2p — > Ti 3d — e g both at the T point. 
The high peak at 5 eV due to O 2p — > Ti 3d-e g at the T point. A third shoulder peak 
at 6 eV is due to the transition O 2p — > Ti 3d-e g also at the T point. As in the case 
of CaTiO^ , the structure in the high frequency part has both lower heights and are 
shifted to higher frequencies. The cause is the same as discussed above. 



BaTi0 3 : 



Lastly, let us look at figures |6|-|8]. If we compare the shape of the imaginary part of the 
dielectric function e%(oS) obtained by our accurate estimate of the transition rate with 
that of we note that agreement with experiment |^2[ is much better when we take 
the energy-frequency dependences of the transition matrix into account. The relative 
weights of the low frequency peaks, at 3.8 eV and 5 eV are correctly reproduced here. 



t i i i i i r 




Photon Energy (eV) 



Figure 6. Comparison of calculated and experimental real part ei(w)(top) and 
imaginary part £2(w) (bottom) of dielectric function of CaTiO^ as a function of Photon 
energy: dotted line is experimental (p5|-p6t), continuous line is theoretical. 



The lower peak is attributed to the transition from the O 2p to the Ti 3c?-t 29 band and 
from the O 2p to the Ti 3d-e g band both at the T points. The next higher peak is the 
because of the transition from the O 2p to the Ti 3d-e g band at the T point. Our present 
study also indicates a peak around 10 eV and the features at higher frequencies follow the 
experimental results closely, although the amplitude seems to have been underestimated 
as compared to the low frequency results. In all cases, although the lower frequency 
part is much better reproduced, the high frequency structures in the theoretical result 
is shifted upwards. As before, we argue that this is probably an artifact of the rigid 
shift of the scissors type approach. A energy-frequency dependent self-energy of the 
type given by the GW method should provide the necessary correction. In addition, the 
wrong heights at higher frequencies could also arise from the fact that we have neglected 
transitions to some of the higher energy conduction bands. A more complete, perhaps 
a three panel calculation at higher energies, should correct this. Alternatively, we could 
use the higher order NMTOs [43|, when they are available, since they span a much larger 
energy range. 
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Figure 7. Comparison of calculated and experimental real part ei(w)(top) and 
imaginary part 62 (w) (bottom) of dielectric function of SrTiO^ as a function of Photon 
energy: dotted line is experimental (^2|), continuous line is theoretical. 

4. Conclusion 

We have proposed here a modified expression for the optical conductivity as a 
convolution of an energy-frequency dependent transition matrix and the energy resolved 
joint density of states. The main motivation was to generalize it to disordered systems, 
where the traditional reciprocal space formulation breaks down due to the failure of 
Bloch's Theorem. In order to be confident in our new formulation we have applied it 
here to the three alkaline earth perovskite titanates in their paraelectric phases. The 
results are in reasonable agreement with experimental data. The agreement is as good 
as we can expect from a LDA calculation. This formulation will now form the starting 
point of a two- fold generalization : first, combining with the ASR to random systems 
and as then as a starting point for a many-body GW formulation. 
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Figure 8. Comparison of Calculated and experimental real part ei(w)(top) and 
imaginary part £2(w) (bottom) of dielectric function of BaTiO^ as a function of Photon 
energy: dotted line is experimental (^2|), continuous line is theoretical. 
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